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ABSTRACT 

The ~800 yr-old pulsar J1846— 0258 is a unique transition object between rotation-powered pulsars 
and magnet ars: though behaving like a rotation-powered pulsar most of the time, in 2006 it exhibited 
a distinctly magnetar-like outburst accompanied by a large glitch. Here we present X-ray timing 
observations taken with the Rossi X-ray Timing Explorer over a 2.2-yr period after the X-ray outburst 
and glitch had recovered. We observe that the braking index of the pulsar, previously measured to be 
n = 2.65 ± 0.01, is now n = 2.16 d= 0.13, a decrease of 18 d= 5%. We also note a persistent increase in 
the timing noise relative to the pre-outburst level. Despite the timing changes, a 2009 Chandra X-ray 
Observatory observation shows that the X-ray flux and spectrum of the pulsar and its wind nebula 
are consistent with the quiescent levels observed in 2000. 

Subject headings: pulsars: general — pulsars: individual (PSR J1846— 0258) — supernovae: individual 
(Kes 75 — X-rays: stars) 


1. INTRODUCTION 

PSR J1846— 0258 is the very young (~800 yr) X-ray 
pulsar at the center of the supernova remnant Kes 75 
(Gotthelf et al. 2000). The pulsar has a spin-period of 
326 ms and usually exhibits properties common to those 
of rotation powered-pulsars (RPPs), including powering 
a bright pulsar wind nebula (PWN). This pulsar is no- 
table for having an unusually high spin-down-inferred 
magnetic field (B = 5 x 10 13 G), and is one of the few 
with a measured braking index (n = vvjv 1 = 2.65d=0.01, 
where v is the spin- frequency, and v and v its derivatives; 
Livingstone et al. 2006). Against expectations, this pul- 
sar exhibited distinctly magnetar-like behavior in 2006 
May-July, when it showed several X-ray bursts, an X- 
ray flux increase (Gavriil et al. 2008), a sizable rotational 
glitch with remarkable “overshoot” recovery (Livingstone 
et al. 2010; Kuiper & Hermsen 2009), as well as spectral 
changes (Kumar & Safi-Harb 2008; Gavriil et al. 2008; 
Ng et al. 2008). This pulsar, evidently a RPP/magnetar 
transition object, presents a unique opportunity to ex- 
plore the long-term relationship between magnetic activ- 
ity and neutron star spin-down. 

Bursts of X-rays and variable X-ray flux, as observed 
in most magnetars and PSR J1846— 0258, are proposed 
to originate from small- or large-scale reorganizations 
of the magnetic field (e.g. Thompson et al. 2002). As 
the magnetic field is intimately connected with the tem- 
poral evolution of pulsars, comparing timing behavior 
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before and after an episode of magnetic activity in a 
neutron star could provide important insight into the 
physics of neutron star magnetospheres. None of the 
magnetars have measured braking indices, despite their 
young ages (r c ~ 1 — 100 kyr) 3 . This is because of a 
significant level of low-frequency timing noise (Gavriil 
& Kaspi 2004) and in several cases, large glitches (Dib 
et al. 2008, and references therein), which prevent a mea- 
surement of n in every known magnetar. By contrast, 
PSR J1846— 0258 rotates relatively steadily, allowing for 
a measurement of n with regular monitoring observations 
with the Rossi X-ray Timing Explorer {RXTE). Between 
its discovery in 1999 (Gotthelf et al. 2000) and the out- 
burst and glitch in 2006, it experienced a level of timing 
noise typical of young RPPs along with one small glitch 
(Av/v ~ 2.5 x 10 9 near MJD 52210) and one small 
candidate glitch {Avjv < 5 x 10 -8 near MJD 52910). 

Also available from regular RXTE monitoring is the 
pulsed flux of PSR J1846— 0258. As reported by Gavriil 
et al. (2008), the pulsed flux had returned to its quiescent 
value roughly two months after the 2006 May outburst 
and no further flux variations have been observed (Liv- 
ingstone et al. 2010). However, because RXTE is a non- 
focusing X-ray telescope, no information about the total 
flux or the phase-averaged spectrum are available from 
these data. To measure these quantities and to confirm 
the RXTE results, a focusing X-ray instrument such as 
the Chandra X-ray Observatory is required. 

Indeed, Chandra observations of the pulsar and PWN 

3 See http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 
for a catalog of known magnetars. 
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revealed flux and spectral changes at the time of the 2006 
outburst (Gavriil et al. 2008; Kumar & Safi-Harb 2008; 
Ng et al. 2008). The pulsar’s total flux rose consider- 
ably, but equally as interesting was its change in spec- 
trum. The quiescent spectrum of PSR J1846— 0258 is 
much like that of other young, high spin-down luminos- 
ity RPPs: a simple power law. Chandra observations 
showed that while in outburst, the spectrum softened 
significantly such that it became reminiscent of those ob- 
served from Anomalous X-ray Pulsars (AXPs), namely, 
well described by a power law with an additional ther- 
mal component. In addition, the superior angular res- 
olution of Chandra allowed detailed observations of the 
PWN, showing that a marginal increase in flux may have 
occurred between 2000 and immediately post-outburst 
(Kumar & Safi-Harb 2008). The effect of magnetar-like 
outbursts on PWNe is an open question, as none of the 
bona fide magnetars power nebulae (although there is 
some debate about the possible existence of a PWN sur- 
rounding the AXP IE 1547.0—5408, see Vink & Bomba 
(2009) and Tiengo et al. (2010)). Kumar & Safi-Harb 
(2008) suggested a causal relation between a possible in- 
crease in the PWN flux observed in 2006 and this or 
past magnetar-like outbursts of the pulsar, while Kar- 
galtsev & Pavlov (2008) suggested that the PWN was 
over-luminous compared to those of other young pulsars, 
perhaps owing to previous, unseen magnetar outbursts. 
However, a revised distance estimate of ~6 kpc (Leahy 
& Tian 2008), rather than the previously claimed 19 kpc 
(Becker & Helfand 1984), reduces the implied nebular 
X-ray efficiency to 77 = L PWN /E ~ 0.02. While this is 
still large, it is similar to that observed from the Crab, 
hence need not be powered by previous outbursts. Thus 
any long-term effect of the 2006 outburst on the PWN 
could help clarify this point. 

In this paper, we report on 2.2-yr of RXTE timing 
observations of PSR J1846— 0258 in the post-magnetic 
activity, post-glitch recovery period. We perform phase- 
coherent and partially phase-coherent timing analyses af- 
ter the 2006 glitch had largely recovered, and we report 
a post-burst measurement of n = 2.16 d= 0.13, smaller 
than the pre-outburst value at the 3.8cr level. We also 
quantify the increase of the timing noise observed over 
the bursting episode and discuss the implications of these 
observations. In addition, we report on 2009 Chandra ob- 
servations of the pulsar and associated PWN which show 
that the pulsar and PWN flux and spectra are consistent 
with their pre-outburst values. 

2. RXTE OBSERVATIONS AND ANALYSIS 

RXTE observations of PSR J1846— 0258 were made 
using the Proportional Counter Array (PCA; Jahoda 
et al. 1996; Jahoda et al. 2006). The PCA consists of 
an array of five collimated xenon/methane multi-anode 
proportional counter units (PCUs) operating in the 2 - 
60 keV range, with a total effective area of approximately 
6500 cm 2 and a field of view of ~ 1° FWHM. 

The RXTE data of PSR J1846— 0258 are unevenly 
spaced over llyr from 1999 April 18 through 2010 April 
22 (MJD 51286 - 55308). Observations taken between 
1999 April 18 - 21 (MJD 51286 - 51289) are excluded 
because they cannot be unambiguously phase connected 
with the rest of the data. Data from 2000 January 31 - 
2008 December 10 (MJD 51574 - 54810) were reduced 


and analyzed previously and details can be found in Liv- 
ingstone et al. (2006) and Livingstone et al. (2010). Data 
taken between 2009 January 27 and 2010 April 22 (MJD 
54858 - 55308) are described here for the first time. 

Data were collected in “GoodXenon” mode, which 
records the arrival time (with l-jas resolution) and en- 
ergy (256 channel resolution) of every unrejected event. 
Typically, two to three PCUs were operational during an 
observation. We used the first Xenon layer of each op- 
erational PCU and extracted events in channels 4-48 
(approximately 2 - 20keV), as this produces good qual- 
ity profiles for individual observations. 

Observations were downloaded from the HEASARC 
archive 4 and data from each active PCU were merged 
and binned at (1/1024) s resolution. Photon arrival times 
were converted to barycentric dynamical time (TDB) at 
the solar system barycenter using the J2000 source posi- 
tion RA = 18 h 46 m 24!94 d= 0!01, Decl = -02°58 / 30.1 // ± 
0.2" (Helfand et al. 2003) and the JPL DE200 solar sys- 
tem ephemeris. 

The phase-coherent ephemeris from Livingstone et al. 
(2010) was used to fold the time series for each new obser- 
vation with 16 phase bins. Resulting profiles were cross- 
correlated with a high significance template profile, pro- 
ducing a single Time-Of-Arrival (TOA) for each obser- 
vation. This process produced 265 TOAs between MJD 
51574 and 55308 with a typical uncertainty of ~ 10 ms 
(~3% of the pulse period). The TOAs were fitted to a 
timing model using the pulsar timing software package 
TEMPO 5 . 

3. TIMING ANALYSIS AND RESULTS 
3.1. Phase- coherent Timing Analysis 

In order to make a significant measurement of a deter- 
ministic value of v and thus n, we restricted our phase- 
coherent timing analysis to MJD 54492 - 55308 (2008 
January 27 - 2010 April 22), because earlier observations 
are highly contaminated by glitch recovery and timing 
noise, as discussed in Livingstone et al. (2010). 

We obtained a single phase-coherent timing solution 
(with no phase ambiguities) fitting only v and if for the 
100 TOAs during this time period, shown in the top panel 
of Figure 1. The residuals show a very significant phase 
contribution from a second frequency derivative (i.e. ~15 
phase turns over roughly 2 yrs). We therefore added v 
to our phase-coherent fit, and the resulting residuals are 
shown in the middle panel of the Figure. The fitted v 
corresponds to a braking index of n = 1.888 =b 0.002. 
As visible in the middle panel of Figure 1, significant 
timing noise remains in the data; the timing residuals 
are not Gaussian distributed. As a result, the formal 
la uncertainty on n from this global phase-coherent fit 
underestimates the true uncertainty. Spin parameters 
from this fit are given in Table 1. 

We fitted higher order frequency derivatives to 
“whiten” the phase residuals, a common procedure to 
lessen the contaminating effect of timing noise on fit- 
ted parameters (e.g. Kaspi et al. 1994). Fitting 12 fre- 
quency derivatives (the maximum possible given current 
machine precision) removes the majority of the timing 
noise, though does not render the phase residuals entirely 

4 http://heasarc.gsfc.nasa.gov/docs/axchive.html 

5 http: / / www. atnf . csiro . au / research / pulsar /t empo/ 
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Gaussian, as shown in the bottom panel of Figure 1. Ta- 
ble 2 shows the variation of n as derivatives are fitted, as 
well as x 2 !/* The value of n varies between 1.89 and 2.95 
as higher order derivatives are fitted, without converging 
to a single value, rendering the true value of n ambiguous 
from this analysis. Nevertheless, the range of measured n 
values from this analysis is relatively narrow: the timing 
noise has increased and clearly contaminates, but does 
not completely dominate if. In cases where a parameter 
is dominated by a noise process, it can be several orders 
of magnitude larger and often of the wrong sign (e.g. 
Hobbs et al. 2010), neither of which are seen here. In 
cases where timing noise contaminates a measurement of 
a deterministic parameter but may not dominate, as for 
these PSR J1846— 0258 data, using a partially coherent 
timing analyis can be useful to find the true value. 

3.2. Partially Coherent Timing Analysis 

To mitigate the effects of timing noise, we performed 
a partially coherent timing analysis. We created 48 
short phase-coherent timing solutions spanning 2000- 
2010 (MJDs 51574-55308). For each short timing so- 
lution, we fit only v and v. The time span included in 
each subset was determined from the requirement that 
the reduced x 2 of the fit was ~1, and that no red noise- 
like structure was visible in the data. In addition, we 
allowed the data included in each timing solution to 
overlap by ~l/2, where sufficiently dense sampling was 
available. This improves coverage, which is of particu- 
lar importance while the timing noise level is very high. 
Figure 2 shows the resulting v measurements spanning 
2000 - 2010. This analysis excludes observations taken 
between 2006 May 31 and 2007 January 27, when glitch 
recovery and timing noise prevented a coherent timing 
solution. As discussed in Livingstone et al. (2006) and 
confirmed by our current analysis, from 2000 - 2006 May, 
if increased very regularly, except at the small glitch in 
2001. The steady increase in z> corresponds to a braking 
index of n = 2.65 ± 0.01, as measured from a phase- 
coherent analysis of these data. The large glitch (visible 
as a sudden decrease in if by ~3% in Fig. 2) and the 
increase in timing noise had largely recovered by the be- 
ginning of 2008, as shown in the Figure. 

To obtain a post-burst measurement of n, we ignored 
all timing data prior to MJD 54492 (2008 January 27), 
where the aforementioned glitch recovery and timing 
noise dominate. We performed a weighted least-squares 
fit to 16 v measurements spanning MJDs 54492 - 55308 
(2008 January 27 - 2010 April 22, shown in the in- 
set of Figure 2). Given the large scatter in the post- 
burst if measurements, and the known effects of timing 
noise on these data, it is likely that the formal uncer- 
tainties significantly underestimate the true uncertain- 
ties. Thus, to better estimate the uncertainty on z>, 
we employed a bootstrap error analysis, helpful when 
formal uncertainties may underestimate the true un- 
certainties (Efron 1979), and previously used for this 
same purpose in Livingstone et al. (2005a, b). This re- 
sults in if = 3.13(19) x 10 _21 s“ 3 , corresponding to 
n = 2.16T0.13, where the uncertainty from the bootstrap 
estimate is larger than the formal uncertainty by a factor 
of ~2.4. This new measurement of n is smaller than the 
pre-outburst value of n = 2.65 d= 0.01 at the 3.8a level 
(or 9.1a if only the formal uncertainties are considered). 


Thus, the braking index decreased by An = — 0.49±0.13, 
following the period of magnetar-like activity in 2006. 
This is the first observed significant measurement of a 
change of a braking index. 

To further confirm a change in n, we performed an 
identical partially coherent timing analysis on the pre- 
outburst measurements of if. Because of the small glitch 
and candidate glitch, the data were separated into three 
large subsets. The first subset (prior to the small glitch 
near MDJ 52210), contains only three measurements of 
z>, resulting in n = 2.63 ± 0.04 (where uncertainties for 
this measurement are the formal uncertainties because a 
bootstrap error analysis can not be performed with no 
additional degrees of freedom). The second subset lies 
between the small glitch and a 78-day gap in the data 
(containing a candidate glitch). In this data subset, we 
performed a weighted least squares fit and bootstrap er- 
ror analysis on 7 measurements of if. This resulted in 
n = 2.61 =b 0.07. We repeated this analysis for the 8 
measurements of if prior to the magnetar-like outburst, 
resulting in n = 2.68 =b 0.03. The three values of n are in 
agreement with each other and the value obtained from 
a fully phase-coherent timing analysis (n = 2.65 ± 0.01; 
Livingstone et al. 2006). All pre-outubrst measurements 
of n are systematically larger than that measured post- 
outburst. Note that the uncertainties on the two mea- 
surements of n from partially coherent analyses (with 
bootstrap uncertainties) to pre-outburst data are smaller 
than the uncertainty for the post-outburst value of n, 
despite similar data spans fitted in each case. This is in- 
dicative of an increase in the timing noise post-outburst, 
discussed further below. 

A complicating factor in timing some magnetically ac- 
tive neutron stars is that pulse profile changes often ac- 
company radiative changes and/or glitches (e.g. Kaspi 
et al. 2003). This can affect measured timing parame- 
ters and must therefore be quantified. It has previously 
been reported that no significant changes in the pulse 
profile were detected during the outburst (Livingstone 
et al. 2010; Kuiper & Hermsen 2009). We further ver- 
ified that the pulse profile in the ~2-yr period immedi- 
ately prior to the outburst is not statistically different 
from the summed profile from the ~2-yr of data used to 
measure n, shown in Figure 3. 

In order to further analyse the timing noise contami- 
nating these data, we performed a second partially co- 
herent timing analysis, this time with each data subset 
having z^, z>, and if fitted. The same conditions of x 2 *, ~ 1 
and Gaussian-distributed residuals, were applied to de- 
termine the length of each subset, and once again, the 
period between 2006 May 31 and 2007 January 27 was 
excluded owing to the lack of a coherent timing solution. 
Six values of v were obtained before the outburst, while 
nine values of v were obtained after the outburst, shown 
in Figure 4. The measurements of v before the outburst 
indicate the regular rotation of the pulsar during this pe- 
riod, while the single value of v above the average (visi- 
ble in the inset) occurs directly after the candidate glitch 
near MJD 52910, providing the best evidence for a glitch 
during this period. As visible in the Figure, the value of v 
changed dramatically immediately after the outburst, to 
a maximum of ~200 times the quiescent value, as well as 
varying in sign, indicating a dramatic increase in timing 
noise during the period of glitch recovery. The magnitude 
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of v decays as the glitch recovers during 2007. The inset 
of Figure 4 shows the variation of v on a smaller scale. 
While the effects of glitch recovery and timing noise have 
subsided substantially by 2008, the post-outburst varia- 
tion of if remains larger than its pre-outburst behavior. 

3.3. Timing Noise 

Qualitatively, the timing noise in the 2.2-yr period used 
to obtain the post-burst measurement of n is larger than 
that observed prior to the outburst, though much smaller 
than in the initial aftermath of the outburst, when no 
phase-coherent timing solution was possible. One mea- 
sure of the change in timing noise can be found by com- 
paring the RMS residuals from MJDs 54492 - 55308 
with those in a similar time span from before the out- 
burst. Fitting v, z>, and if phase-coherently to the 2.2- 
yr segment of data spanning MJD 53086 - 53879, just 
before the outburst, results in a timing solution with 
RMS residuals of 11.4ms (0.035 periods), a factor of 
~5.5 smaller than the RMS residuals from MJDs 54492 - 
55308, of 63.6ms (0.19 periods). This shows that the 
timing noise post-burst is significantly larger than be- 
fore the magnetar-like outburst. This is reflected in the 
new measurement of n by the increased uncertainty com- 
pared to that on the pre-outburst value of n for similar 
measurement baselines. 

Another measure of the increase in timing noise is the 
time span that can be included in each partially coherent 
measurement of z> shown Figure 2. When the pulsar is 
less noisy, more data can be included in each short tim- 
ing solution while satisfying the conditions that x 2 ^ ~ 1 
and that no red noise-like structure remain in the data. 
The pulse profile and pulsed flux are steady, important 
because variability in either could cause changes in TOA 
uncertainties, and thus affect the time span for each if 
measurement. From 2000 until 2006 May, on average, 
each measurement of z> was obtained with data spanning 
111 zb 26 days, while in 2007 (when glitch recovery was 
still a significant effect), each measurement of z> spanned 
an average of 33 zb 20 days. From 2008 - 2010, the aver- 
age time span for each measurement was 68 =L 16 days. 
Thus, nearly four years post-outburst, the pulsar remains 
noisier than prior to the outburst. 

A well known measure of timing noise is the A 8 param- 
eter, defined as the contribution to the rotational phase 
of the pulsar from a measurement of if over a period of 
10 8 s assuming that if is entirely dominated by timing 
noise (Arzoumanian et al. 1994). This parameter is of 
limited value for a pulsar where if is dominated instead 
by secular spin-down, where most of the phase contribu- 
tion from if is due to magnetic braking or another deter- 
ministic spin-down mechanism. To quantify the change 
in timing noise observed in PSR J1846— 0258, we define 
a similar parameter which quantifies the contribution to 
the rotational phase from the measurement of the third 
frequency derivative, i>‘, over a time span of ~ 2.5 x 10 7 s. 
The time span is optimized for this particular pulsar, as 
it is the approximate amount of time required to obtain 
a significant measurement of z>‘, while allowing several 
measurements to be made given the available data span. 
Thus, in analogy with the A 8 parameter we define: 

A « 


We measured the A - parameter for each ~ 2.5 x 10 7 s 
segment of data where a phase-coherent timing solution 
was available, and show the results in Figure 5. The value 
of A- increases dramatically after the 2006 outburst, 
after which it decays, but by 2010 has not returned to 
the pre-outburst quiescent level. 

4. CHANDRA OBSERVATIONS AND ANALYSIS 

PSR J1846— 0258 was previously observed with the 
Chandra X-ray Observatory in 2000 and 2006 with 37 ks 
and 155 ks exposures, respectively. Serendipitously, the 
latter observation was carried out 7 days after the out- 
burst. We obtained a new Chandra imaging observa- 
tion with a 44.6 ks ACIS-S exposure (ObsID 10938) on 
MJD 55053 (2009 August 10), over 3 years after the 
outburst. The observation was taken in 1/8 subarray 
mode, which gives a short frame time of 0.4 s, result- 
ing in negligible (< 3%) pileup of the pulsar counts. 
All the data reduction and analysis were performed with 
Cl AO 4.1 6 . Figure 6 shows the exposure-corrected im- 
ages of PSR J1846— 0258 and the associated PWN in 
the 1 - 7keV energy range. While the PWN exhibits 
no obvious change in the overall morphology between 
the three observations, time variabilities are observed in 
some small-scale features. The most prominent area of 
variability is the northern clump located ~ 7" northeast 
of the pulsar. We extracted the count profiles of the 
clump from the exposure-corrected images using 6" wide 
boxes (as indicated in Fig. 6), and show the results for 
each epoch in Figure 7. The clump morphology evolved 
from a single peak in 2000 to a double peak in 2006, and 
back to a fainter single peak in 2009. However, there is 
no evidence for bulk motion of the clump, as is visible in 
the Figure. 

For the spectral analysis, we extracted the spectrum 
of the entire PWN excluding the central 3" radius, and 
fitted it with an absorbed power-law model. The col- 
umn density was fixed at AT H = 4.0 x 10 22 cm -2 dur- 
ing the fit, to provide a direct comparison with previ- 
ous studies (Ng et al. 2008; Kumar & Safi-Harb 2008). 
In the 2009 observation, we found a photon index of 
T = 1.90 zb 0.03 and an absorbed flux of /o^ s _io = 
(1.37 =b 0.05) x 10 -11 ergss" 1 cm -2 in the 0.5 - lOkeV 
range (hereafter, quoted uncertainties are 90% confidence 
levels). This is fully consistent with the PWN spec- 
trum in 2000, which has T = 1.88 d= 0.03 and /o^ s _io = 
(1.39 =b 0.02) x 10 -11 ergss" 1 cm" 2 (Ng et al. 2008). 

For the pulsar itself, counts were extracted from a 
2" radius aperture and placed in spectral bins with 70 
counts each. We followed Ng et al. (2008) to account for 
the nebular contamination using a power-law model with 
fixed T = 1.9 and /o^f-io = 1-0 x 10“ 12 ergss -1 cm -2 . 
The pulsar spectrum in 2009 is adequately fitted by an 
absorbed power-law model with T = 1.1 zb 0.1. Because 
the above fit provides a reasonable description of the 
data, no blackbody component is required. However, 
in order to compare these data with previous studies we 
also fit an absorbed blackbody plus power-law model. We 
find that the thermal component is not statistically sig- 
nificant, but we derive an upper limit on the blackbody 
temperature of 0.25 keV at the 90% confidence level, sig- 
nificantly lower than the value of 0.9zb0.2keV measured 

6 http://cxc.harvard.edu/ciao4. 1/ 
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during outburst. Table 3 compares the best-fit spec- 
tral parameters among the three epochs, indicating that 
the pulsar flux in 2009 is at a similar level as in 2000, 
much lower than that in 2006 during the outburst. For 
completeness, we also fitted a two blackbody model to 
the 2009 Chandra observation. While the goodness-of- 
fit is statistically similar to that obtained when fitting 
a power-law model, the two blackbody model gives an 
unphysically high temperature of 1.6 ± 0.1 keV for the 
hotter thermal component. In addition, there is evidence 
for a continuous power-law component from 1 — 300 keV 
(Kuiper & Hermsen 2009), rendering the two blackbody 
model of limited interest for this pulsar. 

5. DISCUSSION 

5.1. Timing Noise and the Braking Index 

We have observed a decrease in the braking index of 
PSR J1846— 0258 from n = 2.65 ± 0.01 to n = 2.16 ± 
0.13, corresponding to An = —0.49 =t 0.13, or a decrease 
of 18 ± 5%. The change in n was accompanied by an 
increase in the timing noise of the pulsar, which remains 
larger than the pre-outburst level, nearly four years after 
the glitch and outburst on 2006 May 31. Previous long- 
term observations of n in young pulsars have shown that 
they are remarkably steady. In PSR B1509— 58, timing 
observations over 21 years show that n varies by only 
~ 1.5% (Livingstone et al. 2005b), while the Crab pulsar 
exhibits variations on the order of 5% (Lyne et al. 1993). 

There are two possible ways to interpret the measure- 
ment of An = —0.49 d= 0.13. The first is that the true 
n decreased by a significant amount after the 2006 out- 
burst. The second is that the increased timing noise is 
causing an apparent decrease in n. We discuss each of 
these interpretations next. 

5.1.1. Variable Braking Index 

If the true value of the braking index changed perma- 
nently at the time of the magnetar-like outburst, what 
could be the physical cause of such an effect? From the 
spin-down law derived assuming magnetic dipole braking 
(e.g. Ostriker & Gunn 1969), 

—87 r 2 B 2 R 6 sm 2 a 3 /ri . 

u =—z 7 ?— (2) 

we can infer that d 2 I/dt 2 > 0 7 , da/dt > 0 or dB/dt > 0 
will cause n < 3. 

While it is hard to imagine a physical situation caus- 
ing an accelerated growth or decay in 7, varying values 
of a or B have been considered in the past. Counter- 
alignment of the magnetic field (i.e. an increasing a) 
results in n < 3, even if An = 0. A sudden increase in 
the rate of change of a would produce An < 0. How- 
ever, this is difficult to invoke for the observations of 
PSR J1846— 0258 because the pulse profile shows no vari- 
ation over the relevant time period (e.g. Figure 3). A 
small change in a could still be possible if our line of 

7 Note that the different relationship between I and An than 
for B and a arises from correctly including the variability of I in 
the spin-down luminosity, E , when deriving the above spin-down 
law. The full form of the spin-down luminosity is E = AkIvv + 
47 u/ 2 dl/dt, where only the first term is considered if I is constant. 
Using the full form of E in the derivation of the spin-down law 
leads to the dependence of d 2 I/dt 2 on An. 


sight is crossing the center of the pulsar beam, as a large 
change in the pulse profile may not be required from a 
small change in field orientation. However, this is hard to 
reconcile with the lack of detected radio pulsations from 
the source (Archibald et al. 2008), which is typically in- 
terpreted as our line-of-sight missing the magnetic pole, 
from where the radio emission is thought to originate, 
and crossing only the wider X-ray beam. 

An increase in the magnitude of H, without a change 
in the orientation of the field could also cause a decrease 
in n (Blandford & Romani 1988; Blandford et al. 1983). 
Making all the assumptions of perfect dipole spin-down 
but allowing dB/dt > 0, a braking index of n = 2.65 
implies a time scale of growth for the magnetic field of 
~8000yr, while n = 2.16 shortens the growth timescale 
to ~3500yr. The possibility of magnetic field growth is 
intriguing given the magnetic activity observed prior to 
the change in n. Bhattacharya & Soni (2007) suggest 
that magnetar-strength fields emerge over a period of 
time, as shielding currents dissipate. If the effective B is 
currently in such a period of growth, the smaller value of 
n could result (Lyne 2004). 

The standard spin-down law (Eq. 2) is a major idealiza- 
tion even for a rotation-powered pulsar; for magnetars, 
the picture is almost certainly much more complicated. 
According to one version of the magnetar model, there is 
a one-to-one relationship between n and the large-scale 
twist angle between the north and south hemispheres of 
magnetic field, A0 N -s (Thompson et al. 2002). Here, a 
decrease in n would imply an increase in the twist angle. 
A braking index of n = 2.65 corresponds to a twist angle 
of A (j) ~ lrad, while n = 2.16 implies a twist angle of 
A cj) ~ 2 rad. If the magnetic field remains approximately 
constant before and after the outburst, such an increase 
in the twist angle should be accompanied by an increase 
in the X-ray luminosity of ~ 50%, whereas no persistent 
increase in L x is observed in PSR J1846— 0258 as evi- 
denced by the consistent flux value measured with Chan- 
dra in 2000 and 2009. Further, while the above model as- 
sumes a global magnetic field twist, Beloborodov (2009) 
argues that such global self-similar twists do not present 
a viable explanation for magnetar behavior. Instead, he 
suggests that magnetospheric currents are confined to 
narrow regions on the most extended field lines, which 
are not responsible for the bulk of magnetar X-ray emis- 
sion. If this picture is correct, the spin-down of the star 
can vary without accompanying changes in X-ray lumi- 
nosity. In addition this model predicts that an increase 
in spin-down torque (though not necessarily monotonic) 
should occur sometime after a radiative event, and should 
eventually return to the pre-outburst torque value. Qual- 
itatively, this seems to present an explanation of the ob- 
served timing variability in PSR J1846— 0258, however, 
it provides no quantitative prediction for n. 

A variable n provides an unambiguous test of the the- 
ory of spin down presented by Melatos (1997). He posits 
that the radius relevant to neutron star spin down is 
not the point-like neutron star radius, but a somewhat 
larger “vacuum radius” where field aligned flow breaks 
down. This radius is large enough that the system can 
no longer be treated as a point dipole, resulting in mod- 
ifications to the standard spin-down predictions. In the 
context of this model, a measurement of z>, and a 
uniquely predicts n. While there is no estimate of a for 
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PSR J1846— 0258, there are now two measurements of n 
so the theory can be tested. Given the observed change 
in n, and assuming that a is stable over the magnetar- 
like outburst, this theory predicts that the magnetic field 
should have increased by a factor of ~6. This is not 
observed, however, as the magnetic field estimate has 
increased by just 0.3% compared with the pre-outburst 
value, in contradiction to this theory. Alternatively, if a 
were allowed to change at the time of the event, Melatos 
(1997) predicts that the angle between the spin and mag- 
netic axes would have changed from ~9° to ~4°. Such a 
change might be visible as differences in the pulse profile 
for some lines of sight, but cannot be excluded given the 
lack of observed profile changes. 

Another possible explanation put forth for a static 
measurement of n < 3 is that magnetic field lines are 
deformed due to plasma in the magnetosphere (Bland- 
ford & Romani 1988). A sudden increase in the amount 
of plasma in the magnetosphere of PSR J1846— 0258, 
perhaps injected at the time of the magnetar-like out- 
burst, could cause n to decrease. The best evidence 
for a plasma-filled magnetosphere affecting pulsar spin- 
down is found from the “intermittent” radio pulsar, 
PSR B1931+24, which has dramatic, quasi-periodic 
changes in v correlated with radio pulsations which turn 
on and off (Kramer et al. 2006). Harding et al. (1999) 
propose that spin down can arise from a combination of 
magnetic dipole radiation and wind losses. An increase 
in losses from the wind relative to dipole radiation will 
manifest as a smaller value of n. The spin-down formula 
given by Harding et al. (1999) implies a braking index of 


2v Lp^BR? 
v 2lV6c* ’ 


( 3 ) 


where L p is the kinematic luminosity of the wind, which 
can, in turn, be estimated with a measurement of n. For 
PSR J1846— 0258, a change in n from 2.65 to 2.16 corre- 
sponds to nearly an order of magnitude increase in the 
persistent particle luminosity from L p ~lx 10 36 ergss -1 
to L p ~ 6 x 10 36 ergss -1 . Because the outflowing parti- 
cles travel near the speed of light, the additional particles 
would populate the PWN and be observable as a factor of 
~6 increase in the PWN flux in the 2009 Chandra obser- 
vation as compared to the observation in 2000. However, 
no such flux increase is detected (see §4). The lack of a 
PWN flux increase seems to refute the idea that an in- 
crease in wind losses is responsible for An < 0, however, 
a more rigorous derivation of the relationship between 
wind losses and spin down may provide further insight 
into this issue. 

Interestingly, because a change in the plasma condi- 
tions in the magnetosphere might also affect magneto- 
spheric torques, thus could possibly explain an increase 
in timing noise (e.g. Cheng 1987). As shown by Lyne 
et al. (2010), timing noise in some pulsars can be traced 
to magnetospheric fluctuations. For these radio pulsars, 
there is a correlation between torque variations and pulse 
shape changes. This is difficult to apply in the case of 
PSR J1846— 0258 because no radio pulsations have been 
detected (Archibald et al. 2008) and there is no evidence 
for profile variability in the X-ray band (see Fig. 3). 

Contopoulos & Spitkovsky (2006) note that if the co- 
rotation radius of the magnetosphere, r co _ rot , is less than 


the light cylinder radius owing to imperfect reconnec- 
tion within the magnetosphere, 1 < n < 3 will result. 
If the co-rotation radius decreased at the time of the 
magnetar-like outburst from magnetic field lines open- 
ing, n would also decrease. They parametrize the rela- 
tionship between r co - r ot and n as 


^*co-rot — RLC 


(0 


( 4 ) 


If the initial spin-period of PSR J1846-0258 was Po = 
lms, the implied co-rotation radius prior to outburst 
was r co _ rot = 0.3 6Rlc when n = 2.65 ± 0.01, and 
would have decreased to 0.09 Rlc post-outburst when 
n = 2.16 =b 0.13. Changes to the extent of the co-rotation 
region of the magnetosphere might be visible as changes 
in the pulse profile, however, no such changes have been 
observed from PSR J1846— 0258. In the context of this 
model, however, it is impossible to rule out changes to 
the extent of the magnetosphere as no independent esti- 
mate of the initial spin period is available. In the context 
of this model, taking Pq = lms as the lower limit on the 
birth spin period for PSR J1846— 0258, provides a lower 
limit on the co-rotation radius for any measured value of 
n. Furthermore, Po = 1 ms provides an upper limit on 
the change in co-rotation radius at the two epochs, since 
this quantity decreases with slower birth spin periods. 


5.1.2. Timing Noise 

The second possible interpretation of An < 0 is that 
the true n is constant but remains masked by the high 
timing noise. The increase in timing noise could arise as a 
result of changes to the superfluid interior brought on by 
the unusual 2006 glitch or changes in the magnetosphere 
after the outburst. Though a bootstrap error analysis 
was employed in order to better account for the effect 
of the increase in timing noise, a definitive test is not 
possible. Continued timing observations may be able to 
solve this issue, if the timing noise level continues to de- 
crease and the new braking index remains steady. How- 
ever, it is also possible that the increased level of timing 
noise and the decreased braking index are connected, for 
example, via an increase in the magnetospheric plasma 
density. In that case, if the pre-outburst conditions are 
eventually reobtained, both the timing noise and braking 
index should relax to their pre-outburst values, rendering 
a temporary value of n < 2.65 ambiguous in nature. 

Fluctuations in pinned superfluid in the pulsar inte- 
rior is one of the possible causes of timing noise (e.g. 
Alpar et al. 1986). Thus, one possible explanation for 
the increase in timing noise is that significant changes 
were imparted to the neutron star interior at the time of 
the 2006 glitch and outburst. The glitch was followed by 
an unusual over-shoot recovery (Q ~ 8.7 on a timescale 
of Td ~127days) and a permanent increase in the mag- 
nitude of v (with fractional magnitude Aif/if ~ 0.0041; 
Livingstone et al. 2010). This recovery is thus far unique 
and the origin and long-term consequences of such behav- 
ior is not well understood. By contrast, the permanent 
change in if following the glitch is not unusual when com- 
pared to those measured after other glitches. In addition, 
we note that the change in if is not responsible for the 
observed decrease in n. Because the fractional increase 
in v is three orders of magnitude smaller than the frac- 
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tional change in £>, it is the latter effect that dominates 
the change in n. 

In addition to the detected glitch recovery and perma- 
nent increase in z>, non-monotonic variations in v were 
observed in the aftermath of the glitch. While glitch 
recovery (i.e. a temporary increase in v) and discrete 
jumps in v accompanying a glitch are both established 
phenomena, to our knowledge, no other rotation-powered 
pulsar has experienced changes in timing noise similar to 
those observed from PSR J1846— 0258. However, vari- 
able spin-down torque has been observed in several mag- 
netars. The AXP IE 1048.1—5937 twice experienced ap- 
proximately year-long periods of rapid v variations, i.e. 
sudden, but temporary increases in timing noise (Gavriil 
& Kaspi 2004; Dib et al. 2009). Similar variations in 
if were observed in the Transient AXP XTE 1810—197 
after its 2003 outburst (Camilo et al. 2007). Thus the 
observed change in timing noise in PSR J1846— 0258 can 
be interpreted as yet another example of magnetar-like 
properties from this rotation-powered pulsar, even if the 
phenomenon is currently unexplained. 

5 . 2 . Phase-averaged pulsar flux and the pulsar wind 

nebula 

The Chandra results show that the phase-averaged pul- 
sar flux and spectrum in 2009 have returned to the qui- 
escent values observed in 2000. This agrees with the 
pulsed flux history from RXTE , where the initial flux in- 
crease in 2006 was observed to decay exponentially with 
1/e 55 days (Gavriil et al. 2008), and has since re- 

mained at the quiescent level (Livingstone et al. 2010). 
Given that the 2009 Chandra data were taken over 1000 
days after the 2006 outburst, our findings are not unex- 
pected. The time variability of the small-scale features in 
the PWN could be attributed to magnetohydrodynamic 
instabilities in the flow, similar to what has been ob- 
served in the PWNe powered by the Vela pulsar and PSR 
B1509— 58 (Pavlov et al. 2003; DeLaney et al. 2006) and 
is unlikely to be related to a sudden deposition of par- 
ticles at the time of the outburst. Gavriil et al. (2008) 
reported an energy release of ~ 5 x 10 41 (d/6 kpc) 2 ergs 
(2 - 60keV) in the 2006 outburst. In comparison, the to- 
tal energy released in the giant flare from the magnetar 
SGR 1806—20 in 2004 is 2 x 10 46 ergs (Palmer et al. 2005), 
of which 4 x 10 43 ergs went into particle energy (Gaensler 
et al. 2005). The resulting ratio between the particle en- 
ergy and electromagnetic radiation is 2 x 10 -3 . If the 
same ratio holds for PSR J1846— 0258, then the energy 
of the injected particles would be only ~ 10 39 ergs. With 
a 5-field strength of 15/iG in the PWN (Djannati-Ata’i 
et al. 2007), the synchrotron cooling timescale of a par- 
ticle emitting at 5keV is ~300yr. These particles would 
have induced an X-ray flux enhancement in the PWN of 
2 x 10 _17 ergss _1 cm -2 . This is 5 orders of magnitude 
lower than the X-ray flux of the northern clump, and 
6 orders lower than that of the entire PWN, too small 
to be detected. Moreover, if the particles travel isotropi- 
cally with a typical post-shock flow speed of c/3, the flux 
would have already spread over the entire PWN of 10" 
radius. Therefore, we do not expect the 2006 outburst 
to have any observable effects on the PWN. 

6. CONCLUSIONS 


The observed change in n after the magnetar-like out- 
burst in PSR J1846— 0258, if shown to be steady via 
ongoing timing observations, has important implications 
for the physics of neutron star spin-down. A decrease in 
n could have several origins, and a detailed theoretical 
framework is necessary for interpreting this observation. 

Most theoretical descriptions of a changing n require 
an accompanying persistent change in radiative behav- 
ior of the pulsar, while we observe neither pulse profile 
variability or persistent flux enhancement. An increase in 
particle wind losses relative to dipole losses does not pro- 
vide a good description of An < 0 for PSR J1846-0258 
because of the lack of a persistent increase in PWN lu- 
minosity (Harding et al. 1999). However, variability in 
magnetospheric plasma remains a promising avenue for 
future consideration, especially considering the recent 
report of variable spin-down rate correlated with radio 
pulse shape changes for several pulsars (Lyne et al. 2010). 
While no variability in the X-ray pulse profile is detected 
in PSR J1846— 0258, short time scale variability would 
not be detectable in RXTE observations which are typ- 
ically 1.5 to 2 hr long. 

The timing noise in PSR J1846— 0258 is observed to 
be of a higher level than prior to the outburst. That 
is, even four years after the glitch and magnetic activ- 
ity, the pulsar is rotating less regularly than in its pre- 
outburst quiescent state. It is interesting to note, how- 
ever, that the current timing noise in PSR J1846— 0258 
would not be unusual if observed in any young pulsar. 
Rather, it is the sudden change in the level of timing 
noise in PSR J1846— 0258 that is noteworthy. Since the 
timing noise is simply of a higher level and not other- 
wise different from that observed in other pulsars, the 
phenomenon cannot be used as a diagnostic of previous, 
unseen magnetic activity in other pulsars. 

The observed decrease in n and increase in timing noise 
reported here may or may not be permanent. Regu- 
lar monitoring observations beyond the RXTE era may 
help to answer this question, as well as to search for 
future magnetar-like X-ray outbursts and glitches from 
PSR J1846— 0258. 

The 2009 Chandra observations of PSR J1846— 0258 
show that the total flux and spectrum are consistent with 
the quiescent values observed in 2000. No significant 
variability was detected in the PWN that can be asso- 
ciated with the 2006 outburst, nor, given the energetics 
of the outburst, would any be expected. The variability 
of the PWN observed in the 2006 Chandra observation 
is most likely unrelated to the outburst, and is probably 
similar in origin to the variation in small-scale features 
seen in other PWNe. 
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Fig. 1. — Timing residuals of PSR J1846— 0258 spanning MJDs 54492 - 55308 (2008 January 27 - 2010 April 22). The top panel shows 
residuals with v and if fitted. The middle panel shows residuals with z/, V and i> fitted while the bottom panel shows residuals with 12 
frequency derivatives fitted. 
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Fig. 2. — Evolution of the frequency derivative of PSR J1846— 0258 over ~10yrs of RXTE timing observations. Measurements of i > 
overlap by ~l/2 where sufficient data are available. The effect of the two confirmed glitches are visible on the Figure, the first near 
MJD 52210 as a small discrete jump in if, the second as the large decrease in if near MJD 53886 (note that the full effect of the glitch on 
if is not shown here in the interest of making visible the smaller changes in if at other epochs, see Livingstone et al. (2010) for additional 
details about the large glitch). The inset shows measurements of if from MJDs 54492 - 55308, with the best-fit slope shown as a solid line 
and the =bl<j uncertainties (from the bootstrap analysis which better accounts for the increase in timing noise) shown as dotted lines. The 
slope corresponding to the pre-outburst n is shown as a dashed line (colored blue in the online version), where uncertainties are roughly 
an order of magnitude smaller than those from the post-outburst era, so are not visible on the Figure. 
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Fig. 3. — Pulse profile of PSR J1846— 0258 for two sections of data of ~2yr. Two cycles are shown for clarity. The solid line shows 
the pulse profile before the outburst in 2006 May. The dotted (red in the online version) line shows the summed pulse profile for 2008 
January - 2010 April. Subtracting the two profiles results in residuals with x 2 ^ ~ 1.3, where the probability of this x 2 *, or higher occurring 
by chance is 19%. Thus the two profiles are not statistically significantly different. 
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Fig. 4. — Second derivative measurements for PSR J 1846— 0258 from 2000 - 2010. The three measurements occurring after the glitch 
and outburst are 1-2 orders of magnitude larger than the subsequent measurements and vary in sign, indicating that they are severely 
contaminated by glitch recovery and/or timing noise. The inset shows measurements of i> on a smaller scale to highlight the smaller 
variation away from the glitch recovery. The variation in i> during the period from 2008 to 2010 is larger than in the pre-outburst era, 
and the mean value is systematically smaller. The one value of i> pre- burst that is significantly larger than the average is directly after the 
candidate glitch near MJD 52910, and is the best evidence that a glitch actually occurred at that epoch. 
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Fig. 5. — A quantification of the timing noise in PSR J1846— 0258 over 10 yr. Each point is a measurement of the A — parameter for 
approximately 2.5xl0 7 s. This provides an estimate of the amount of timing noise observed in the pulsar, and shows a dramatic increase 
after the large glitch observed in 2006. 


2000 

2006 

2009 







v JKL\ » . 



Ijfrairfc 




• 

I ■ 

\ \ 



\ \ 

\ r d \ 4 


N 

10" 




0 1.3e-7 3. 9e-7 9.0e-7 1.9e-6 4.0e-6 8.1e-6 1.6e-5 3.3e-5 6.5e-5 1.3e-4 


Fig. 6. — Exposure-corrected Chandra ACIS images of the PWN associated with PSR J1846— 0258 in 1 - 7keV. The 6" wide box along 
a position angle of 27° (North through East) is analyzed in detail for each epoch in Figure 7. 
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Fig. 7. — Zoom-in of Figure 6 showing the northern clump of the PWN (boxed region in Figure 6) in the 1 - 7keV range. The bottom 
right panel illustrates the count profiles extracted from the 6"- wide box (as shown in Figure 6), indicating the evolution of the clump 
between epochs. 





Post-Outburst Observations of PSR J1846— 0258 


15 


TABLE 1 

2008 - 2010 Spin Parameters for PSR J1846-0258. 


Phase-coherent timing analysis 

Date range (Modified Julian Day) 

54492.089 - 55308.598 

Date range 

2008 Jan 27 - 2010 Apr 22 

Number of TOAs 

100 

Epoch (Year) 

2009 Mar 1 

Epoch (MJD) 

54834.0 

i/ (Hz) 

3.0621185502(4) 

z> (10 11 s 2 ) 

-6.664350(2) 

i> (10- 21 s- 3 ) 

2.725(3) 

Number of derivatives fitted 

2 

RMS residuals (ms) 

63.6 

Partially phase- coherent timing analysis 

i> (10“ 21 s -a ) 

3.13(19) 

Braking index (n) 

2.16(13) 


Note. — Quoted uncertainties are the formal 1<t uncertainties as 
reported by TEMPO for the fully phase- coherent timing analysis. 
Details about the uncertainty on i> and n for the partially coherent 
timing solution are from a bootstrap analysis, as explained in the 
text (see § 3.1). 


TABLE 2 

Variation of n with number 

OF FITTED FREQUENCY 
DERIVATIVES 


Derivatives 

fitted 

n 

xl 

2 

1.888(2) 

39.71 

3 

2.010(3) 

6.12 

4 

1.980(8) 

6.00 

5 

2.05(1) 

5.48 

6 

2.08(2) 

5.50 

7 

2.51(4) 

3.51 

8 

2.60(4) 

3.26 

9 

2.95(7) 

2.93 

10 

2.91(7) 

2.73 

11 

2.1(1) 

2.12 

12 

2.1(1) 

2.14 


Note. — Braking index varia- 
tion with number of fitted deriva- 
tives. The number of degrees of 
freedom for two fitted derivatives 
is 96. Uncertainties are the for- 
mal 1(7 uncertainties returned by 
TEMPO for the timing solution 
spanning MJD 54492 - 55308. 
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TABLE 3 

SPECTRAL FITS TO PSR J1846-0258 


Epoch 

JV H 

(10 22 cm- 2 ) 

r 

j?abs,PL 
J 0.5—10 
(10 -12 ergs 
cm -2 s _1 ) 

Minabs,PL 
J 0.5-10 
(10 -12 ergs 
cm -2 s -1 ) 

fcT 

(keV) 

f abs,BB 
J 0.5—10 
(10 -12 ergs 
cm -2 s -1 ) 

punabs,BB 
■10.5-10 
(10 -12 ergs 
cm -2 s -1 ) 

xS/dof 

2000 

2006 

2009 

4.0* 

4.0* 

4.0t 

i n+0-8 
iU -0.3 

1.9 ±0.1 
1.0 ±0.1 

4.2 ±0.2 

1.3 ±0.3 
3.6 ±0.5 

6.1 ±0.4 
31 ±6 

5.2 ±0.7 

0.4 ±0.4 
0.9 ±0.2 
0.16 ±0.09 

o.oi±8;8? 

1.7 ±0.2 
0.010 ± 0.005 

o.i±8;? 

3.2 ±0.4 
0.30 ±0.15 

0.94/38 

0.99/134 

0.84/72 


References. — The 2000 and 2006 results are from Ng et al. (2008). 
Note. — The uncertainties quoted are 90% confidence intervals, 
t - held fixed in the fit. 



